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ABSTRACT 


The  optimum  gain  o£  uniformly  spaced  arrays  of  isotropic  sources 
or  dipoles  is  investigated  theoretically  in  this  paper.  The  formulation 
is  processed  with  the  aid  of  an  array  matrix.  The  optimum  gain  and  the 
corresponding  excitation  are  expressed  directly  in  terms  of  the  elements 
of  the  array  matrix. 


1522-1 


iii 


Page  1  of  37 


ON  THE  OPTIMUM  GAIN  OF  UNIFORMLY  SPACED 
ARRAYS  OF  ISOTROPIC  SOURCES  OR  DIPOLES 


INTRODUCTION 

For  a  uniformly  spaced  array  consisting  of  a  given  number  of 
isotropic  sources,  the  renowned  synthesis  of  Dolph,1  including  the  ex¬ 
tension  by  Riblet,  2  offers  the  minimum  beamwidth  for  a  prescribed 
side-lobe  level.  It  is  not  an  optimum  design  from  the  point  of  view  of 
maximizing  the  directivity  or  the  gain.  Ma  and  Cheng3  have  recently 
presented  another  synthesis  which  optimizes  the  gain  under  the  con¬ 
dition  of  a  prescribed  side-lobe  level.  Because  of  the  polynomial  formu¬ 
lation,  neither  of  these  two  methods  can  conveniently  be  applied  to  arrays 
of  directive  sources,  such  as  those  consisting  of  dipoles.  The  problem 
that  deals  solely  with  the  maximum  or  the  optimum  gain  of  uniformly 
spaced  arrays  was  first  investigated  by  Uzkov4  in  a  very  elegent  formu¬ 
lation.  By  means  of  an  orthogonal  transformation  in  vector  space,  he 
obtained  some  important  results  concerning  the  optimum  gain  of  an 
array  of  isotropic  sources.  In  particular,  he  showed  that  the  optimum 
gain  of  an  end-fire  array  of  isotropic  sources  as  the  separation  approaches 
zero  is  numerically  equal  to  N2  ,  where  N  denotes  the  number  of  sources. 
He  also  showed  that  when  the  separation  is  equal  to  X  /2,  the  optimum 
gain  is  numerically  equal  to  N.  Although  he  indicated  that  the  method 
can  be  applied  to  arrays  of  directive  sources,  he  did  not  elaborate.  The 
excitation  of  the  arrays  to  produce  the  optimum  gain  was  not  discussed 
in  his  work.  Several  years  later.  Block,  Medhurst  and  Pool5  proposed 
an  optimization  method  based  upon  the  impedance  matrix  defined  for  the 
elements  of  an  array.  Only  a  very  limited  calculation  was  reported  in 
their  work.  Recently,  Stearn6  made  some  extensive  calculations  on  the 
gain  of  an  end-fire  array  of  half-wave  dipoles,  and  its  corresponding 
excitation  for  various  separations  based  upon  their  method. 

In  this  work,  we  shall  formulate  the  problem  with  the  aid  of  an 
array  matrix.  The  method,  in  certain  respects,  is  equivalent  to  the 
impedance  method  except  that  all  the  essential  results,  namely,  the 
gain  and  the  corresponding  excitations,  can  be  expressed  directly  in 
terms  of  the  elements  of  an  array  matrix.  The  formulation  is  particularly 
effective  in  dealing  with  arrays  of  short  dipoles  as  well  as  isotropic 
sources.  In  the  present  report,  we  shall  consider  only  broadside  arrays, 
and  ordinary  end-fire  arrays.  The  latter  are  characterized  by  a  restraint 
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that  the  progressive  phase  shift  between  adjacent  elements  is  assumed 
to  be  equal  to  the  electrical  distance  between  them.  This  type  of  end-fire 
array  does  not  yield  the  optimum  gain,  because  of  the  restraint  on  the 
phase  shift.  We  treat  these  arrays  here  because  the  formulation  is  very 
much  like  that  for  broadside  arrays.  In  a  subsequent  report,  a  general 
treatment  of  end-fire  arrays  will  be  given  without  imposing  the  above- 
mentioned  restraint  on  the  progressive  phase  shift. 


FORMULATION 


As  an  illustration  of  the  method,  the  simplest  case  of  a  broadside 
array  consisting  of  2n  isotropic  sources  as  shown  in  Fig.  1(a)  will  be 
treated  here.  The  field  pattern  of  such  an  array  is  given  by 


(1) 


F{9) 


'  I  Ai  c0S[(‘-i)D 

i=l 


cos  6 


where 

D  =  2rrd/X 

d  =  separation  between  adjacent  elements 

0  =  angle  of  inclination  measured  from  the  axis 

of  the  array 


A^  =  amplitude  of  the  i-th  pair  of  elements 
2n  =  N  =  number  of  elements . 


The  amplitude  A^  is  assumed  to  be  real,  but  it  may  be  either  positive 
or  negative.  The  corresponding  power  pattern  can  then  be  written  in 
the  form 


(2) 


S(0)  =  [F(0)]2  = 


l 
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The  directivity  or  the  gain  of  the  array,  designed  to  be  broadside,  is 
determined  by 


i=l  j=l 


where 


-  TT 

■ 

'  /  V 

/  \ 

(4) 

Pij  Pji  =  2 

r  cos 

'  o 

(  i  -  JD  cos  9 

cos 

fj  -  -  Djcos  0 

=  -  [si.j(D)  +  Si+j.jtD)]. 
2 


The  sr(D)  function  appearing  in  (4)  is  defined  by 


(5) 


sr(D)  = 


_  sin  rD 


rD 


r  t  0; 


and 


s0(D)  =  1  , 


r  =  0. 


To  optimize  the  gain,  we  set 


(6) 


p  =  1,  2, 


n 
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which  yields  the  following  set  of  equations: 


(7) 


or 

(8) 


i  i*  a 

i=l  j=l  i=l 

p  =  1,  2,  • •  •  n; 


1 


j=l 


n. 


Since  the  quantity  on  the  right  side  of  (8)  is  independent  of  p,  we  may 
denote  it  by  K,  and  (8)  can  be  written,  in  the  matrix  form, 


(9)  [p]  {A}  =  {K}  . 

The  elements  of  the  square  matrix  [p],  which  is  symmetric,  are  given 
by  (4).  For  convenience,  [p]  will  be  designated  as  the  array  matrix. 
Denoting  the  adjoint  of  [p]  by  [b],  one  has 


(10)  {A}  =  [p]"1  {K}  =  |  p  f 1  [B]  {K} 


where  j  (3  |  denotes  the  determinant  of  the  array  matrix.  From  (10), 
we  conclude  that  for  optimum  gain,  the  amplitude  distribution  must 
satisfy  the  following  relations: 


n  n  n 

(11)  A,  :  Aj  :  ...  A„  =  7  Blj:  y  Bjj  :  •••  Y  Bnj  . 

j=l  j=l  j=l 
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In  view  of  (8),  Eq.  (?)  may  be  written  as: 


(12) 


PpjAj; 


p  =  1 ,  or  2,  •  •  •  or  n 


where  the  capital  letter  Gj^  signifies  the  optimum  value  of  gj^j.  By 
making  use  of  (11),  and  a  multiplication  rule  in  determinant  theory, 

(12)  can  be  converted  into  the  following  expression: 

n  n 

(13)  GN=  l  y  B../|p|  . 

i=l  j=i 

Equations  (11)  and  (13)  constitute  the  principal  result  of  this  formulation. 
The  simplicity  of  these  formulae  is  manifest  since  only  the  elements  of 
the  array  matrix  are  contained  in  tnese  expressions.  For  the  case  of  a 
broadside  array  with  an  odd  number  of  elements  equal  to  2n-l,  Fig.  1(b), 
Eqs .  (11)  and  (13),  are  still  valid,  except  that  the  elements  of  the 
array  matrix  are  changed  to 


(14)  Pij  ~  j  (si_j  +  si+j_2)  ,  for  N  =  2n  -  1. 

The  method  described  here  can  be  applied  equally  well  to  end-fire  arrays, 
or  arrays  of  directive  sources.  In  the  following  section  we  list  the  ex¬ 
pression  for  the  power  pattern  and  the  corresponding  array  matrix  for 
several  types  of  arrays  composed  of  isotropic  sources  or  of  short  dipoles. 
It  is  assumed  that  the  individual  short  dipoles  have  a  figure-eight  pattern 
described  by  sin  9^,  where  9^  denotes  the  polar  angle  measured  from 
the  axis  of  the  dipole.  For  completeness,  the  cases  discussed  previously 
are  also  included. 
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POWER  PATTERN  AND  THE  CORRESPONDING  ARRAY 
MATRIX  FOR  VARIOUS  TYPES  OF  ARRAYS 


where 

s0(D)  =  1,  sr(D)  =  £HL£P  ,  D  =  —  d. 

rD  X. 


b)  N  =  2n-l 


An  A2  A(  Aj  A2  A 

(a) 


An  A2  2A|  A2  An 

(b) 


Fig.  I.  Broadside  arrays  of  isotropic  sources  (a)  N  =  2n,  (b)  N  =  2n-1. 

The  designation  for  the  elements  of  an  array  with  odd  number 
of  elements  also  applies  to  other  types  of  arrays  sketched  in 
the  subsequent  figures. 
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In  the  following  tabulation,  the  functions  Se(9)  and  S0(9)  have  the  same 
meaning  as  those  defined  for  broadside  arrays  of  isotropic  sources. 

(II)  Broadside  array  of  parallel  dipoles  (Fig.  2) 

a)  N  =  2n 


where 


S(9,cj>)  =  (1  -  sin2  9  cos2  <$>)  Se(9) 
Pij=j  [pi-j(D)  +Pi+j-l<D>] 


Po(D)  =  |  , 


Pr(°) 


sin  rD  cos  rD 

Id"  + 


b)  N  =  2n  -  1 


S(9,4»)  =  (1  -  sin2  9  cos2  4>)  S0(9) 
Pij  =  i  tPi.j(D)  +Pi+j.2(D)]- 


Fig.  2.  A  broadside  array  of  parallel  dipoles. 
The  angles  9  and  £  appearing  in  the 
text  correspond  to  the  ones  shown  in 
this  figure . 
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(Ill)  Broadside  array  of  collinear  dipoles  (Fig.  3) 


a)  N  =  2n 


where 


S(0)  =  sin2  0  Se(0) 

Pij  =  ii-j(D)  +  *i+j-l(D> 


f0<D)  =  j,  ir(D) 


sin  rD  cos  rD 
r3  D3  r2  D2 


b)  N  =  2n  -  1 


S(0)  =  sin2  0  So(0) 

Pij  =  +  ^i+j-2^)* 


Fig.  3.  A  broadside  array  of  collinear  dipoles. 


(IV)  Broadside  array  of  crossed  dipoles  (Fig.  4) 

In  this  case  the  two  crossed  dipoles  are  assumed  to  be  of  equal 
amplitude  but  excited  90°  out  of  phase.  The  field  due  to  a  single  unit  in 
the  broadside  direction  (0  =  7r/2,  4>  =  tt/2)  is  therefore  circularly  polari¬ 
zed. 


a)  N  =  2n 

S(0,4>)  =  —  (1  +  sin2  c(>  sin2  0)  Se(0) 

Pij  =  5  [Ci_j(D)  +  Cj+j.jID)] 
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where 


C0(D)=|,  Cr(D) 


0 +  -M  - 


in  rD  cos  rD 
rD  "  r2  D2 


b)  N  =  2n  -  1 

S(9,  <j>)  =  L  (1  +  sin2  4>  sin2  9)SQ(9) 
Pij  -  i  [Ci.j(D)  +Ci+j.2(D)J. 


A2  Aj  A|  A2 


An 


Fig.  4.  A  broadside  array  of  crossed  dipoles; 
the  excitation  of  the  first  pair  is: 

A1X  =  Aj  and  Alz  =  jAj  ,  etc.  The 
angles  9  and  <j>  are  shown  in  Fig.  2. 


(V)  End-fire  array  of  isotropic  sources  (Fig.  5) 

As  explained  in  the  introduction,  the  end-fire  arrays  considered 
here  are  of  the  ordinary  type,  with  a  progressive  phase  shift  between 
adjacent  elements  numerically  equal  to  D. 


a)  N  =  2n 


Se'(9) 

Pij 


n 

■  )  A;  COS 

\r  \ 

[i  -  i_)D(l  -  cos  9) 

> 

IV  2/  \ 

li=l 

L  x  ' 

% 

\  ts2(i-j)(D)  +  s2(i+j-l)(D)] 


where  the  sr(D)  function  is  the  same  as  that  defined  in  Case  (I).  We 
may  remark  that  S2r(D)  can  be  decomposed  as  follows: 
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=  cos  rD  sr(D). 


s2r<D)  =  cos  rD 


sin  rD 
rD 


b)  N  =  2n-l 


Sc(9)  =  i 


V  2 

cos  [(i-l)D(l-cos  0)] 


i=l 


Pij  "  1  [s2(i-j)(D)  +  S2(i+j-2)(D)J  * 


Fig.  5.  An  end-fire  array  of  isotropic 
sources,  D  =  2ir(d/\  ). 


(VI)  End-fire  array  of  parallel  dipole  (Fig.  6) 


a)  N  =  2n 

S(0,  <j>)  -  (1  -  sin2  0  cos2  4>)  Sg^O) 

Pij  =  |  [li-j<D)  + 

where  Se#(0)  is  the  same  as  that  defined  in  Case  (V-a)  and 

1  \  sin  rD  +  cos  rD 

r2  D2  J  rD  r2  D2 


qr(D) 


cos  rD 


(‘ 


=  cos  rD  pr(D) . 
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The  function  pr(D)  appears  in  Case  (II). 
b)  N  =  2n-l 

S( 8 , 4> )  =  (1  -  sin2  9  cos2  cj>) 


Fig.  6.  An  end-fire  array  of  parallel  dipoles; 

the  angles  9  and  <j>  are  shown  in  Fig.  2. 

(VII)  End-fire  array  of  crossed  dipoles  (Fig.  7) 

The  crossed  dipoles  are  assumed  to  be  of  equal  amplitude  but 
excited  90°  out  of  phase. 

a)  N  =  2n 

S(9)  =  —  (1  +  cos2  9)  Sje) 

2  e 

Pij  =  \  K-j  +  qi+j-l]  ' 

Although  the  pattern  for  this  case  is  different  from  that  of  an  end-fire 
array  of  parallel  dipoles,  the  array  matrix,  and  hence  the  optimum 
gain,  is  the  same  for  both  cases. 
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b)  N  =  2n-l 


s(0)  =  I{i  +  cos2  0)  s0'(e) 

2 

Pij  =  \  Ili-j  +  1i+j-2l  • 


-Jf — — -J& 


U- 


^ 


An6 


.  2n-l  _ 
j-r-D 


A,e 


3D 
1  2 


A,e 


<§ 


A,e 


JT 


A9e 


•3D 

j  2 


Ane 


Fig.  7. 


An  end-fire  array  of  crossed  dipoles;  the  excitation 
of  the  first  unit  a1j_the  right  is: 

~j  2 


.D 
J  2 


MX 


=  Aj  e  ,  A^  =  jAj  e  ,  etc. 


The  angle  0  and  4>  are  shown  in  Fig.  2. 


NUMERICAL  COMPUTATION 

The  numerical  computations  for  the  amplitude  distribution  and  the 
optimum  gain  based  upon  (11)  and  (13)  seem  to  be  straightforward. 

This  is  indeed  the  case  if  the  order  of  the  array  matrix  is  equal  to  two, 
or  if  D,  the  separation  between  the  elements,  is  not  too  small.  When 
n  is  greater  than  two,  and  D  is  less  than  tt  for  the  broadside  arrays  or 
less  than  tt  / 2  for  the  ordinary  end-fire  arrays,  the  computation  becomes 
rather  difficult  even  with  the  aid  of  an  IBM-7  090  computer.  The  fact 
that  such  a  giant  machine  could  not  evaluate  accurately  a  3  x  3  determi¬ 
nant  or  invert  the  associated  matrix  without  going  into  multiple  precision 
programming  is  truly  unexpected.  It  should  be  mentioned  that  D  =  tt 
for  the  broadside  cases,  and  D  =  tr/2  for  the  ordinary  end-fire  cases, 
correspond  to  the  demarcations  of  the  Dolph  synthesis  and  the  Riblet's 
extension.  For  convenience,  we  shall  call  the  region  of  D  lying  below 
these  two  values  the  "super-gain"  region.  Pending  a  more  accurate 
evaluation,  we  shall  present  here  only  the  data  of  Gj^  for  the  values  of 
D  where  the  accuracy  was  certain.  These  results  are  plotted  in 
Fig.  8-13.  In  regard  to  the  corresponding  amplitude  distributions,  the 


1522-1 


12 


data  are  too  numerous.  Figures  14-19  give  a  few  sample  curves  for  some 
typical  distributions.  In  the  case  of  broadside  arrays  of  isotropic  sources, 
as  discussed  in  the  following  section,  it  is  possible  to  find  the  limiting 
values  of  Gjyj  as  d  -*•0.  The  dotted  lines  of  the  Gj^  curves  in  the  super-gain 
region  correspond  to  the  interpolated  values  between  these  limiting  values 
and  the  accurately  computed  values.  For  other  cases,  where  these 
limiting  values  have  not  yet  been  found,  only  the  accurately  computed 
values  are  plotted.  Each  set  of  curves  contain  also  a  plot  of  Gmax/G00 
and  where  Gmax  and  Gmin  denote,  respectively,  the  first  maxi¬ 

mum  and  the  first  minimum  of  the  Gjvj  curves.  G^  denotes  the  asymptotic 
value  of  Gjsj  as  d  —  oo  .  Numerically,  G^  is  equal  to  N  for  arrays  of  iso¬ 
tropic  sources,  and  is  equal  to  3/2  N  for  arrays  of  dipoles.  In  the  case 
of  broadside  arrays  of  isotropic  sources,  we  also  plotted  a  curve. 

Fig.  8(e),  labelled  G0/ G^  where  G0  denotes  the  limiting  value  of  Gj^  as 
d  -*-0.  A  discussion  on  the  values  of  GQ  is  given  in  the  next  section. 

THE  LIMITING  VALUE  OF  Gn  FOR  BROADSIDE 
ARRAYS  OF  ISOTROPIC  SOURCES 

For  N  =  3  or  4,  it  is  relatively  simple  to  evaluate  the  limiting  value 
of  Gn  as  D  approaches  zero.  This  can  be  done  by  expanding  the  relevant 
function  sn{D)  in  a  series  of  D2  and  retaining  the  leading  terms  of 

!  P  I  and  ^  ^  Bij  • 

The  result  gives 


(15) 


lim  G3 
D-0 


lim  G* 
D-0 


2.25. 


When  N  =  5  or  6,  the  same  procedure  calls  for  seven  terms  in  the  series 
expansion  of  sn(D)  which  already  involves  a  very  tedious  algebraic 
manipulation.  The  result  yields 


(16)  lim  G5  =  lim  G6  =(-r  -j  )  =  3.515625  . 

D-0  D-0  \2  4/ 

To  obtain  this  simple  result  we  must  deal  with  many  clumsy  numbers 
such  as  (6)l4/15!  .  It  should  be  mentioned  here  that  these  calculations 
were  performed  before  we  had  access  to  Uzkov's  article.  In  view 
of  the  orderly  figures  contained  in  (15)  and  (16),  we  then  proposed  a 
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conjecture  that  for  arrays  with  odd  number  of  elements  N  or  those  with 
even  number  of  elements  N  +  1  the  limiting  value  must  be  given  by 


(17) 


lim 
D— 0 


l 

6 


After  reading  Uzkov's  paper,  we  realized  that  as  far  as  these  limiting 
values  are  concerned  they  can  most  conveniently  be  obtained  from  his 
method.  Although  Uzkov  gave  these  values  only  for  end-fire  arrays,  it 
can  be  shown  by  means  of  the  orthogonal  transformation  which  he  out¬ 
lined  that,  in  general, 

N-l 

(18)  lim  Gn(0o)  =  7  (2n+l)  (cos  0Q) 

D-0 


where  Pn(cos  0O)  denotes  the  Legendre  polynomial  of  order  n,  and  0O 
denotes  the  angle,  measured  with  respect  to  the  axis  of  the  array,  in 
which  direction  the  array  was  designed  for  the  optimum  gain.  The  values 
of 


lim  Gn(0o) 

D— 0 

are  plotted,  on  db  scale,  in  Fig.  20(a)  and  (b).  For  end-fire  arrays, 
without  the  phase  restraint  mentioned  previously  in  the  introduction, 
one  has 


(19) 


lim  Gj^(0)  = 
D-0 


N-l 

7  (2n+l)  =  N2  . 

n=0 


This  formula  was  originally  enunciated  by  Uzkov.  For  the  broadside 
array,  (18)  reduces  to 


(20) 


lim  Gj^j 
D-0 


N-l 

£  (2n+l)P*(l). 
n=0 


1522-1 


14 


The  same  formula  applies  to  the  ordinary  end-fire  arrays  as  D  -*  0  . 

It  can  be  shown  that  the  sum  given  by  (20)  is  indeed  identical  to  the  one 
given  in  product  form  by  (17).  Looking  back  at  the  hard  work  which  we 
did  in  arriving  at  (17),  one  has  to  admire  Uzkov's  ingenious  method  of 
successive  orthogonal  transformations  in  providing  the  answer  for  these 
limiting  values  of  Gjvj.  The  only  consolation  for  our  innocent  hard  work 
is  that  without  it,  the  alternative  expression  for  (20)  as  given  by  (17) 
probably  would  not  be  recognized  at  first  glance.  In  regard  to  the 
numerical  values  of  the  limiting  values  of  G^(tt/2)  as  given  by  (17)  or 
(20),  it  is  of  some  interest  to  point  out  that  these  can  be  approximated 
quite  accurately  by  an  even  simpler  formula.  From  the  Peirce-Foster 
Table,7  one  finds  that 


(21) 


4n  >  /  3  5  #  _  2n-l  I  >  2(2n-l) 

tt  1  2  4  2n-2y  tt 


hence  the  mean  value  of  4 n/ir  and  (4n-2/rr)  may  be  used  an  an  approxi¬ 
mate  value  of  Gn(7t/2),  i.e.. 


(22) 


._  2(2n-|)  _  2N+1 

TT  TT 


The  expression  given  in  (22)  is  certainly  the  asymptotic  value  of 
Gj^(tt/2)  for  large  values  of  N.  The  fact  that  it  may  be  used  for  all 
values  of  N  in  approximating  (17)  is  clear  from  the  following  table: 


TABLE  Of(-  •  -  •••  — )  AND  ^±1 
\  2  4  N-l  *  tt 


N 

f  3 

5  ' 

N  f 

2N  +  1 

\  2 

4 

2N-1  / 

TT 

3 

2.250 

2.228 

5 

3.515 

3.501 

7 

4.785 

4.775 

9 

6.056 

6.048 

11 

7.328 

7.321 

13 

8.600 

8.594 

15 

9.873 

9.868 

17 

11.14 

11.14 

19 

12.42 

12.42 
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Finally,  we  wish  to  emphasize  that  the  result  found  here  also  agrees  with 
the  conclusion  drawn  by  Bou.wkamp  and  Bruijn8  that  for  a  continuously 
distributed  line  source,  N  -+a>,  the  optimum  gain  is  without  bound.  How¬ 
ever,  for  discrete  arrays  there  is  a  unique  solution. 


COMPARISON  OF  GN  WITH  THE  GAIN  OF  A 
UNIFORMLY  EXCITED  ARRAY 

Before  we  conclude  this  report,  it  is  important  to  point  out  that 
the  gain  of  a  uniformly  excited  array  with  the  same  number  of  elements 
and  spacing,  which  will  be  denoted  by  gj^>  is,  in  general,  slightly  less 
than  the  optimum  gain  except  in  the  super-gain  region.  Figures  21-25 
show  the  comparison  between  g^j1'  and  G^  for  some  typical  cases.  In 
view  of  the  practical  value  of  uniformly  excited  arrays,  we  have  compiled 
a  fairly  complete  set  of  curves  of  g^j1'  for  various  types  of  arrays.  The 
data  will  soon  be  published  as  a  separate  report. 
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Fig.  13.  Optimum  gam 


Fig.  14.  Amplitude  distribution  of  a  broadside  array 
of  six  isotropic  sources. 
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Fig.  15.  Amplitude  distribution  of  a  broadside  array 
of  eight  parallel  dipoles. 
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Fig.  17.  Amplitude  distribution  of  a  broadside  array 
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Fig.  18.  Amplitude  distribution  of  an  ordinary  end -fire 
array  of  four  isotropic  sources. 
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ig.  19.  Amplitude  distribution  of  an  ordinary  end-fire 
array  of  six  parallel  dipoles. 
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Fig.  22.  Comparison  of  g8  and  Gg  of  broadside 
arrays  of  parallel  dipoles. 
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